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An isospin analysis of B — > tttt decays yields sin 2a, where a is the usual CKM angle a = 
aTg[—VtdVf*i,/ (VudV*),)] without hadronic uncertainty if isospin is a perfect symmetry. Yet isospin 
symmetry is broken not only by electroweak effects but also by the u and d quark mass difference 
— the latter drives -k^ — rj,ri' mixing and converts the isospin-perfect triangle relation between the 
B — > TTTT amplitudes to a quadrilateral. The error incurred in sin 2a through the neglect of the 
Qs . resulting isospin- violating effects can be significant, particularly if sin 2a is small. 



I. INTRODUCTION 



In the standard model, CP violation is characterized by a single phase in the Cabibbo-Kobayashi-Maskawa (CKM) 
matrix, rendering its elements complex. The CKM matrix is also unitary, so that determining whether or not this is 
empirically so is a central test of the standard model's veracity lit] . Ascertaining whether the angles of the unitarity 
'^ • triangle, a, (3, and 7, empirically sum to tt and whether its angles are compatible with the measured lengths of its 
\^ ' sides lie at the heart of these tests of the standard model. 

OA ' We study the impact of isospin violation on the extraction of sin 2a from an isospin analysis in B — )■ tttt decays [^ . 

^S^ Isospin is broken not only by electroweak effects but also by the u and d quark mass difference. The latter drives 
^\ Ti^ — 77,77' mixing [y|, which, in turn, generates an amplitude in i? — > tttt not included in the isospin analysis. Thus, 
^^ , although the effect of electroweak penguins is estimated to be small [Q-g , when all the effects of isospin violation are 
rl^ ' included, the error in the extracted value of sin 2a can be significant [Q. 
"^ To review the isospin analysis in i? ^ tttt decays, due to Gronau and London B, let us consider the time-dependent 

Q-i, asymmetry A(t) ||]: 



•''''^tTlSS"»'^"''>-(lTSp)^'°'^'"''' <" 



, ^ where r^, = (V,;V,j/F,iV',i)(>lfc,//lfc,) = e-^'*- A,^JA,^,, A,^, s A(B° ^ fcp), and Am s i3„ - Bi |. 
. . 1 Denoting the amplitudes 5+ — > 7r+7r°, S" -^ tt^tt", and i?° -^ vr+vr^ by A^^ , yl"°, and A'^ , respectively, and 

introducing Aj to denote an amplitude of final-state isospin /, we have 

^-A+- =A2-Ao ■ AO" = 2A2 + Ao ; -i=A+o = iA^ , (2) 

where analogous relations exist for A^", A , and A in terms of A2 and Aq- If isospin were perfect, then the Bose 
symmetry of the J = tttt state would permit amplitudes merely of / = 0, 2, so that the amplitudes B^ -^ tt^tt^ 
would be purely 1 — 2. In this limit the penguin contributions are strictly of A/ = 1/2 character, so that they cannot 
contribute to the 1 = 2 amplitude: no CP violation is possible in the tt^tt^ final states. The penguin contribution in 
B^ —f TT+TT^, or in B^ — > tt+tt^, can then be isolated and removed by determining the relative magnitude and phase 
of the / = to / = 2 amplitudes. We have 
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" "" {A2 - ^0) (1 - ^) ' ^ ' 

where z{z) = ^0/^2(^0/^2) and A2/A2 = exp(-2z0t) with (j)t = arg(T4dK*b) and (/)„i + 0t = /3 + 7 = tt - a in the 
standard model 0]. Given 1^4"' |, |j4°°|, |A+°|, and their charge conjugates, the measurement of Imr^+^- determines 
sin 2a, modulo discrete ambiguities in arg((l — z)/{l — z)), which correspond geometrically to the orientation of the 
"triangle" of amplitudes associated with Eq. (0) , namely 

A+- + 2A°" = V2A+° , (4) 

with respect to |A+°| = |^~"| and that of its charge conjugate. The triangles' relative orientation can be resolved via 
a measurement oilmrj^a^o as well g, and thus sin 2a is determined uniquely. 

II. ISOSPIN VIOLATION AND tt^-ti, rj' MIXING 

We examine the manner in which isospin- violating effects impact the extraction of sin 2a, for isospin is merely an 
approximate symmetry. The charge difference between the u and d quarks engenders a A/ = 3/2 electroweak penguin 
contribution, which is outside the scope of the dehneated isospin analysis Q, although methods have been suggested to 
include them 1 10 1^ . This is the only manner in which the u-d charge difference enters our analysis, so that we term this 



source of isospin breaking an "electroweak effect." The u-d quark mass difference can also engender a A/ ~ 3/2 strong 
penguin contribution through isospin-breaking in the hadronic matrix elements. Moreover, strong-interaction isospin 
violation drives 7r° — ri,ri' mixing [BJ, admitting an / = 1 amplitude. Although electroweak penguin contributions are 
estimated to be small Mfl, other isospin- violating effects, such as iT^-r],r]' mixing, can be important 0,^3]. 

To include the effects of vr"-?], rj' mixing, we write the pion mass eigenstate JTr*^) in terms of the SU(3)/-perfect 
states \(f>3), \(f>s), and |0o): where, in the quark model, j^a) = \uu — dd)/V2, \(j)s) — \uu + dd — 2ss)/^/6, and 
\(j)o) = \uu + dd + ss)/^/3. Explicit relations between the physical and SU(3)/-perfect states can be realized by 
expanding QCD to leading order in l/Nc, momenta, and quark masses to yield a low-energy, effective Lagrangian 
in which the pseudoscalar meson octet and singlet states are treated on the same footing |§,|l3|- Diagonalizing the 
quadratic terms in (/)3, (/)8, and 0o of the resulting effective Lagrangian determines the mass eigenstates tt^, rj, and 77' 
and yields, to leading order in isospin violation H, 



k°) 



|(/)3)+e(cos6i|08)-sin6l|0o>)+e'(sin6i|08)-|-cos6'|0o>) , (5) 



where cos 6*108) — iim0\(f>o) = \ri) + 0(e), and sin6'|08) + cos6'|0o) = W) + 0{e'). Moreover, e — eoXCOsS and 
e' = -2eoXsin6', with ^ = 1 + (4m|. - 3m^ - m^)/(m^ -mj) « 1.23, x = 1/x, eo = V3(md - m„)/(4(ms -m)), and 
TO = (to„ + niii)/2 H. Thus the magnitude of isospin breaking is controlled by the SU(3)-breaking parameter nis — rh. 
The 77-77' mixing angle 9 is found to be sin 26* = -(4V2/3)(to|- - 77^^)/(to^, - to^) so that 6 « -22° ||]. The resulting 
e — 1.14eo is comparable to the one-loop-order chiral perturbation theory result of £ = 1.23eo in 77 ^ •7r+7r~7r'' |l^,p[. 
Empirical constraints also exist on the sign of tt^-t;, 77' mixing. That is, the ratio of the reduced matrix elements in 
K12 decays, namely, K'^ —* 7r°e+i^e and K^ —> 7r^e+z^e, is given by [ p5[ 

cxpt 

= 1.029 ±0.010. (6) 

Using the Lagrangian of Rcf. g] and the quark masses mq{ii = 1 GcV) of Ref. IQ yields 

I •' + 




jf^" 



= l + \/3e8» 1-018 ±0.010, (7) 



where eg is the 03 — <f>s mixing angle, eg = e cos 9 + e' sin 6. Note, for comparison, that the one-loop-order chiral 
perturbation theory result is 1.022 ||l^. In regard to the sin 2a results to follow, it is worth noting that the isospin- 
violating parameters we have adopted appear conservative with respect to the existing experimental constraints. Using 
■niqlp = 2.5 GcV) of Ah et al. ||l^, we find e = 1.4 • 10^^ and e' = 7.7 • 10^'^; we use these values in the subsequent 
calculations. 
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FIG. 1. Reduced amplitudes in B -^ tttt in the factorization approximation with [Nc, k /rrif,] for a) [2,0.5], b) [3, 0.5] (solid 

line) and [3, 0.3] (dashed line), and c) [cxj, 0.5]. Note that Ar = 2A°° /{{GF/V2)iVubV;a), Ar' = A^' /{{GF/V2)iVubV;d), 
and A^" = y/2A~'^/{{GF/V2)iVubV*^). The charged modes are separated into tree and penguin contributions, so that 



At 



T+,^3 + Pn+Tro and A^ 



T _ 4- P 



where P^±^o includes the isosran- violating tree contribution in A^ as 



well. The shortest side in each polygon is the vector defined by the RHS of Eq. ([10|); it is non-zero only in the presence of 
tt"-?;, rj' mixing. 



In the presence of Tr'^-Ty, rj' mixing, the B — > tttt amplitudes become 



^-" = (vr- 



/off I 50 



£8(7r 08 



|H^«|i3-)+eo(7r-0o|H«''|i?" 



A = (0303|H°"|S"') +e8(0308|H^"|i3") +£o(030o|H^"|i?") , 



(8) 
(9) 



where Eq is the 03 — 0o mixing angle, Eq = e' cos9 — esin6'. Note that either of the 7r° mesons in the B^ -^ tt^tt^ 
amplitude can suffer n^-r], rj' mixing; the factor of two associated with this appears as 2 A in Eq. (px|). The B ^ tttt 
amplitudes satisfy 



2A 



00 



/effl nO\ 






(10) 



and thus the triangle relation of Eq. ( H) becomes a quadrilateral. We ignore the relatively unimportant mass differences 
77ijr± — tIttO and rre^i — 771^0 ■ 



III. RESULTS 



We proceed by computing the individual amplitudes using the A_B — 1 effective Hamiltonian resulting from the 
operator product expansion in QCD in next-to- leading logarithmic (NLL) order [|l6| , using the factorization approxi- 
mation for the hadronic matrix elements. In this context, we can then apply the isospin analysis delineated above to 
infer sin 2a and thus estimate its theoretical systematic error, incurred through the neglect of isospin violating effects. 
The effective Hamiltonian W^ for h -^ dqq decay can be parametrized as |16 



K 



cff 



V2 



( '° 

\i=3 



(11) 



where Oi and Og are as per Ref. |Iq|; we also adopt their Wilson coefficients Ci and Cg, computed in the naive 
dimensional regularization scheme at a renormalization scale of ^ = 2.5 GeV |16[. In NLL order, the Wilson coefficients 
are scheme-dependent; yet, after computing the hadronic matrix elements to one-loop-order, the matrix elements of 
the effective Hamiltonian are still scheme- independent |lq|. This can be explicitly realized through the replacement 
{dqq\l-L''^\h) = {Gp / \/^){dqq\[VubV*^{CfO\ + CfO^) - VtbV*^ Y^Zz CfOiWhY""", where "tree" denotes a tree-level 
matrix element and the C°'^ are from Ref. Il6|. The Cf^ are complex |l9| and depend on both the CKM matrix 
parameters and fc^, where k is the momentum transferred to the qq pair in 6 — > dqq decay. Noting Ref. po| we 
use p = 0.12, rj = 0.34, and A = 0.2205 1lqj2l|] unless otherwise stated. One expects m^/4<fc^<m^/2 [|2[; we use 
k^/ml — 0.3,0.5 in what follows. We define the decay constants {tt^ {p)\d^fj,j5u\0) = — i/ttP^ and (0i(p)|?l7^75M|O) = 
—ify^jPfj., and use the flavor content of the SU(3)/-perfect states to relate /? to /^r. Finally, using the quark equations 
of motion with PC AC and introducing a^ = Cf + Cffi/N^ for i odd and Oi = Cf + Cf^i/Nc for i even, the 



B- 



TT (/)3 matrix element in the factorization approximation with use of the Fierz relations is 



TABLE I. Strong phases and inferred values of sin 2a |2| from amplitudes in the factorization approximation with Nc and 



k^r. 



■ 0.5. The strong phase 25truo is the opening angle between the An and A_r+ amplitudes in Fig. 1, whereas 25gl is 



the strong phase associated with the closest matching sin2Q values, denoted (sin2a)GL, from Im r^+^-/Im r^o^o, respectively. 
The bounds |25gqi| and |2fc Qn| o n 2(5truc from Eqs. (2.12) and (2.15) of Ref. |24] are also shown. All angles are in degrees. 
We input a) sin2a = 0.0432 [@||], b) sin 2a = -0.233 (p = 0.2, 77 = 0.35) |||, and c) sin 2a = 0.959 (p = -0.12,77 = 0.34). 



case 


Nc 


2dtruc 


|2<5gqi| 


2(5gqii 


\2Sgl\ 


(sin2a)GL 


a 


2 


24.4 


24.5 


12.5 


13.4 


-0.145/0.153 


a 


3 


24.2 


16.1 


15.3 


15.4 


-0.107/0.133 


a 


00 


23.8 


55.8 


19.4 


17.7 


-0.0595/0.101 


b 


2 


19.6 


22.0 


8.3 


9.4 


-0.399/-0.0213 


b 


3 


19.4 


12.9 


12.3 


12.4 


-0.349/-0.0599 


b 


00 


19.2 


56.2 


17.3 


15.7 


-0.287/-0.104 


c 


2 


28.3 


34.4 


14.8 


4.9 


0.769/0.701 (*,t) 


c 


3 


28.0 


22.8 


17.4 


4.0 


0.662/0.692 (*,t) 


c 


00 


27.6 


42.7 


21.2 


19.5 


0.912/0.967 



* The matching procedure fails to choose a sin 2a which is as close to the input value as possible. 
t The discrete ambiguity in the strong phase is resolved wrongly. 



V2' 



{ 7T-<j>3\H'^\B-) = ^[VubV:,{tUFB-^,{ml-)a,+tflFB-Amlo)a-2) - Vt^V;^ 



X iiUFB-^^iml-)ia4 + aio + 



{mu +md){mb - niu)' 



^IIFb- Amlo) 



(12) 



^^o (oe 



3, , 1 

X (04 + ^^(a? - Og) - -aio -1 -, r- 

2 2 md[mb - md) 



■))] 



)^cf^nO)/(l- 



gV^o+): where we use i^o^^''(0) = 
= Fb-k/V^, and i^s^Q = Fbtt/V^- 



The transition form factors are given by FB-Ti{q ) = {itt-b- ~ '' 
0.33 and Mo+ = 5.73 GeV as per Refs. [^^^. Also Fb03 = Fb„/\/2, Fb^,^ 
Note that the 04 and ae terms, which are associated with the strong penguin operators, only contribute to the 
{t:~ (j)z\'H'^^\B~) matrix element if niu 7^ md or if /? -Fb-ttI'™^!)) 7^ fTrFB-^.^{m^_) — we neglect this latter contri- 
bution as we set mT^± — m.^o. The strong-penguin contributions, which are isospin- violating, explicitly realize the 
induced A/ = 3/2 effect discussed previously, for the amplitude {iT~(j)3\?i'^^\B~), in concert with the amplitudes 
(7r-7r+|H'=*^|BO) and (0303|H'=^|BO), satisfy the triangle relation of Eq. (|). Thus the {Tr^(t)-i\H''^\B-) amplitude can 
still be deemed purely 1 = 2 even if rriu 7^ md, as in Eq. (^^. In the presence of tt'^-tj,?]' mixing, however, the 
B^ -^ n^TT^ amplitude can no longer be purely / = 2, as the RHS of Eq. (10) is non-zero and Eq. (|j) is no longer 
satisfied. 

Numerical results for the reduced amphtudes Ar and An, where Aj^ = 2A / {{G p / V2)iVubV*^) , Aj^ = 



-r + - 



A'' /{{GF/V2)iVubV:^), and A^° = V2A-y{{GF/V2)iVubV:^), with N, = 2, 3, 00 and P/ml = 0.3,0.5 are 
shown in Fig.|l|. A'j^^ and A^^ are broken into tree and penguin contributions, so that A^° = T^+^g + Ptt 



and A, 



T 

J- IT 



±0 



P^-T^a. Note that "P^±^o" is defined to include the isospin-violating tree contribution in A^ as 



well. The shortest side in each polygon is the vector defined by the RHS of Eq. (10). The values of sin 2a extracted 
from the computed amplitudes with Nc — note that N^. is regarded as an effective parameter in this context — and 
k^ /ml = 0.5 are shown in Table |; the results for k'^ /ml = 0.3 are similar and have been omitted. For reference, the 
ratio of penguin to tree amplitudes in p- -^ tt^ttO is \P\/\T\ - (2.2 - '2..7)%\VtbV*i\/\VubV*d\ for A^^ = 2, 3 and fc^ as 
above. Were electroweak penguins the only source of isospin violation, then |P|/|T| '-^ (1.4 — 1.5)%|V(bV'j^|/|Vubyj^|, 
commensurate with the estimate of 1.6% in Ref. [H. 

In the presence of tt^-t;, ry' mixing, the Aj:^ , Aj^ , and Aj^ amplitudes obey a quadrilateral relation as per Eq. (M). 
Consequently, the values of sin 2a extracted from Im rv+7r- and Im r.,^a^a measurements can not only differ markedly 
from the value of sin 2a input but also need not match. The incurred error in sin 2a increases as the value to be 
extracted decreases; the structure of Eq. (||) suggests this, for as sin 2a decreases, the quantity Im ((1 — z)/(l — z)) 
becomes more important to determining the extracted value. It is useful to constrast the impact of the various 
isospin-violating effects. The presence of A/ = 3/2 penguin contributions, be they from rn^ 7^ nid or electroweak 
effects, shift the extracted value of sin 2a from its input value, yet the "matching" of the sin 2a values from the Im 
r^+^- and Im rjrOTrO determinations is unaffected. This arises as the amplitudes in question still satisfy the triangle 
relations implied by Eq. (^j. The mismatch troubles seen in Table | are driven by Tr^-yy, 77' mixing, though the latter 
shifts the values of sin 2a extracted from Im 



' 7r+7i 



as well. Picking the closest matching values of sin 2a in the two 



final states also picks the solutions closest to the input value; the exceptions are noted in Table |. The matching 



TABLE II. Errors in a were Im r^+^- taken to be sin 2a (j(5Q|raw) and from applying the Gronau-London analysis [g] in the 
presence of isospin-violating corrections (|5a|GL) for amplitudes computed in the factorization approximation with Nc = 2 and 
fc^/mf = 0.5. All angles are in degrees. Cases a), b), and c) are defined as in Table I. 



iVc 



Im r. 



|<5a|, 



-0.346 
-0.514 
0.642 



11.3 

8.8 
16.8 



5.4 

5.0 

11.6 (* 



* The discrete ambiguity in the strong phase is resolved wrongly in this case — see Table I. 



procedure can also yield the wrong strong phase; in case c) of Table | with N^ = 2,3, the triangles of the chosen 
solutions "point" in the same direction, whereas they actually point oppositely. If \A^'^\ and \A \ are small [g| the 
complete isospin analysis may not be possible, so that we also examine the utility of the bounds recently proposed by 
Grossman and Quinn pj| on the strong phase 25tmc = arg((l — ■z)/(l — z)) of Eq. (S). The bounds 2(5gqi and 2(5gqii 
given by their Eqs. (2.12) and (2.15) W%, respectively, follow from Eq. (0), and thus can be broken by isospin- violating 
effects. As shown in Table B, the bounds typically are broken, and their efficacy does not improve as the value of 
sin 2a to be reconstructed grows large. 

To conclude, we have considered the role of isospin violation in B ^ tttt decays and have found the effects to be 
significant. Most particularly, the utility of the isospin analysis in determining sin 2a strongly depends on the value to 
be reconstructed. The error in sin 2a from a Im r^+^- measurement grows markedly larger as sin 2a grows small — 
this is the region of sin 2a currently favored, albeit weakly, by phenomenology ^6 21 2^j2^ . The effects found arise in 



part because the penguin contribution in B ^ tt+tt", e.g., is itself small; we estimate |P|/|r| < 9%|V(by(^|/|14f,V^i^^|. 
Relative to this scale, the impact of Tr^-ry, 77' mixing is significant. This is displayed in another way in Table |^. The 
"penguin pollution" in i? — > tt+tt" is such that were no isospin analysis applied, the error in a would be of the order 
of 10° — 20°. The isospin- violating effects in B ^ tt+tt^ suggest that the error in a is still of the order of 5° after 
the Gronau-London g] analysis is applied. Yet, were the penguin contributions in B -^ tttt larger, pressing the need 
for the corrections of the isospin analysis, the isospin-violating effects considered would still be germane, for not only 
would the A/ — 3/2 penguin contributions likely be larger, but the B —^ nrj and B — > nrj' contributions could also be 
larger as well [ p7[ . To conclude, we have shown that the presence of 7r'^-ry,ry' mixing breaks the triangle relationship, 
Eq. (Q), usually assumed [g[ and can mask the true value of sin 2a. 
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